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The evolution of the Schrodinger-cat states in a dissipative parametric amplifier 
is examined. The main tool in the analysis is the normally ordered characteristic 



function. Squeezing, photon-number distribution and reduced factorial moments are 



■ discussed for the single- and compound-mode cases. Also the single-mode Wigner 

t> : 

CN ■ function is demonstrated. In addition to the decoherence resulting from the inter- 

action with the environment (damped case) there are two sources which can cause 
such decoherence in the system even if it is completely isolated: these are the decay 
of the pump and the relative phases of the initial cat states. Furthermore, for the 
damped case there are two regimes, which are underdamped and overdamped. In the 

X 

first (second) regime the signal mode or the idler mode "collapses" to a statistical 
mixture (thermal field). 
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I. INTRODUCTION 



The linear superposition principle is at the heart of quantum mechanics since, using it, 
one can control the properties of the single states making them more or less pronounced. 
The most significant example reflecting the power of such a principle are the Schrodinger- 
cat states jj, which exhibit various nonclassical effects such as squeezing, sub-Poissonian 
statistics and oscillations inphoton- number distribution [2-4 1, even if the original states are 
close to the classical ones J5]. Such type of states (of a jth mode for convenience) can be 
represented as 



\ a )<t>i = N A\ a j) + ex P(*0j)l - 

where |aj) is a coherent state of the jth mode with complex amplitudes otj, <pj 
phase and Nj is the normalization constant having the form 



(1) 

is a relative 



N 



1 



. (2) 

3 2[l + exp(-2|a j | 2 )cos0 j ] v ' 

Specifically there are three values of (fij, namely, 0, ir and 7r/2, for which (Op) reduces to even 
coherent (ECS), odd coherent (OCS) and Yurke-Stoler (YSS) states, respectively. These 
three states will be frequently used in this paper. 

There are several proposals for generating states (pQ) in various nonlinear processes, e.g. 
see {4, 6-10]. Further, the evolution of cat states in several quantum systems has been 



intensively studied 



111 



181 ]. For convenience of the present work we refer to the interac- 
tion of these states with environment llJl-J. [rsl]. where the common goal of such studies is 



the comprehension and description of the decoherence processes (decoherence is the rapid 
transformation of a pure linear superposition state into the corresponding statistical mix- 
ture state) in the system under observation, and the energy loss due to the interaction 
with the reservoir. Most of these papers adopt the Heisenberg-Langevin approach or the 
Markov master equation approach. Although the non-Markovian dynamics is important in 
quantum physics (theory), a unified and compact treatment of the non-Markovian reservoirs 



19], where also the decay and revival of the Schrodinger-cat states are 



dates as late from 
characterized. 

On the other hand, parametric amplifier (nondegenerate) takes a considerable interest in 
quantum optics since it can perfectly generate two-mode squeezing. Recently, this device 



3 



has been supported with the fast progress of new nonlinear crystals and improved laser 



sources, especially femtosecond lasers 



the interference experiments 21 



201 ] and it has been employed in experiments, e.g. in 



the statistical properties of the 



24- 



261 ] for compound-mode case 



22]. An investigation o 
parametric amplifier with losses 23j] and without losses 
have been performed using different techniques when the modes are initially prepared in the 
coherent states. The anticorrelation in this model is an interesting effect 26], where also 
under certain conditions the variance of the photon number is less than the average pho- 
ton number and the photocounting distribution becomes narrower than the corresponding 
Poisson distribution for a coherent state with the same mean photon number. Moreover, 
the sum photon-number distribution can exhibit collapses and revivals 27] (and references 
therein) similar to those typical in the Jaynes-Cummings model (JCM) 28(|; the former is 
in the photon number domain rather than the time domain. So the question we would 
like to address in this paper is the following: what would be the influence of the interfer- 
ence in phase space on these phenomena? In other words, we study the evolution of the 
cat states (CQ) through the parametric amplifier and take into account the interaction with 
environment. The system alone may be described by tracing a more complete description 
over the environment. Actually, the interaction considered here is quite different from that 
considered in [ll|- 17] because there are two operations controlling its behaviour, which are 



the interference in phase space and the entanglement (correlation) between the signal and 
idler modes as well as between the signal-idler system and the environmental oscillators 
(reservoirs) to which the system may be connected. Such entanglements lead to the increase 
of the marginal entropies of the subsystems (as is known, if the density operator describes 
pure states, then the entropy (S) is zero, otherwise S ^ 0) |29|). For the present system 
we treat the single- and compound-mode cases. The main tools in our treatment are the 



norma, 
is 



30 



ly ordered characteristic functions. The Hamiltonian which governs the interaction 



3l| 



H = ujjOjCij — hg[aia 2 exp(iuit 

3=1 



+ H.C.1 



(3) 



2 oo 



j=ll=X 

j = 1,2, are the annihilation (creation) operators assigned to the signal 
and idler modes, respectively; Uj are the natural frequencies of oscillations of the uncoupled 



where ci,-, (aj), 
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modes with oj = cu\ + U2,oo being a pump frequency (under the assumption of resonance 
frequency); <fi is the initial phase of the pump and H.c. means the Hermitian conjugate term 
to the previous one; g is the real coupling constant including the amplitude of the pump 
(gain coefficient); bji and b^ are the boson annihilation and creation operators of the reservoir 
oscillators, respectively, with the frequencies <pji, and kji are the coupling constants of the 
system-reservoirs interaction. Hamiltonian fl3]) represents three wave mixing in a nonlinear 
crystal described by the second-order susceptibility ■ Specifically, a pump beam with 
frequency u travelling through the crystal creates photons of the signal and idler modes with 
different frequencies cuj, j = 1,2, such that uj = u>i + uj 2 - The whole system is considered to 
be interacting with its surroundings (e.g., via imperfect cavity mirrors), which are modelled 
as reservoirs. 

For completeness, the well-known solution of the Heisenberg-Langevin equations for the 
Hamiltonian d3J) is 30|, 131 ] 



A x {t) = A(*)fli(o) + / 2 (*)4(o) + E[^(o)r 1/ (t) + bUo)r' u (t)}, 

1 

(4) 

00 

A 2 (t) = f 3 (t)a 2 (0) + h(t)a\(0) + £[&L(0)r 2 ,(t) + b 2l (0)T' 2l (t)l 

1 

where Aj(t) = aj(t) exp(iujt); %(0) and 6^(0), j = 1,2, are the initial operators of 
the modes and reservoirs, respectively. The explicit forms for the dynamical coefficients 
fj(t), Tji(t) and r^(t) are given in the appendix A. Further, it is worth mentioning that the 
Langevin forces have the forms 



= ~i ^2 % & ^( ) exp(-iLpjit), j = 1, 2. 



(5) 



These forces satisfy the following commutation rule 



[L J (t),L],(t')]= l3 5 lj/ 5(t-t')l 



(6) 



where 7^ is the cavity decay rate of the jth mode (for more details about the properties of 



the reservoir oscillators, see 



25 



32|). Finally it is worth mentioning that including lossy 



mechanism, the difference mean-photon number between the signal and idler modes in the 
system becomes nonconservative. 
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This paper is organized as follows: In section 2 the basic equations and relations, such as 
two-mode normally ordered characteristic function, quadrature squeezing, Wigner function, 
reduced factorial moments and photon-number distribution are given when the modes are 
initially prepared in the cat states (CQ). In sections 3 and 4 discussion of the results of the 
single-mode and compound-mode cases are performed, respectively, and the conclusions are 
summarized in section 5. 



II. BASIC RELATIONS AND EQUATIONS 

In this section we give the main relations which will be used to investigate the properties of 
the system under consideration. Actually, the calculations are lengthy, but straightforward, 
and for this reason, we give briefly general calculations in such a way that particular results 
can be obtained using suitable choice of parameters. We start by discussing the initial density 
matrix of the system. We assume that the system and reservoirs are initially independent 
and noninteracting before switching on the interaction, and the interaction between them 
starts at t = 0. This means that the density operator could be written initially as a direct 
product [33J: 

p(0)=p f (0)®p r (0), (7) 

where p r (0) is the density matrix of the system of damping oscillators, which we assume 
to have flat (i.e. constant function of frequency) and broad-band reservoir spectra, so that 
the mean number of reservoir quanta (phonons) in the mode I is (bji(0)bji(0)) = (njd) 
independently of I, where the subscript d denotes broad-band reservoir. Also we assume that 
the reservoirs form chaotic systems with mean photon numbers (njd) = [exp(^^-) — 
where the reservoirs-oscillators are at the temperature T, with the frequency (pj, Kb is 
Boltzmann's constant and h Planck's constant divided by 2ir. We proceed assuming that 
Pf(0) is the field density operator, which in our case, where the modes are assumed to be 
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initially in superposition states (OQ), has the form 

Pf(0) = la}^ |o;)^ 2 ^ 2 (a|0 1 (o;| 



(8) 



= N?Ni[p M (0) + psi(O) + pai(O)]. 
The concise form in the second line of ([8]) means that the values of ocj , <pj , and <j>^ should 
be specified to obtain the exact form of the density operator in the first line. Here /3m (0) 
denotes the statistical mixture part of the density operator as 

Pm(0) = |ai)|a 2 )(a 2 |(a 1 | + |ai)| - a 2 )(-a 2 |(ai| 

(9) 

+ | — ai)|a; 2 )(a 2 |(— a%\ + | — an)| — a 2 )(— a 2 \(— «i|; 
Psi(0) denotes symmetric interference part of the density operator in which the two modes 
are in off-diagonal basis of coherent states and can be written as 

Psi(0) = |exp[i(0x + 2 )]| - - a 2 )(a 2 |(ai| 

(10) 

+ exp[— + 2 )]|q;i}|q! 2 }(— a 2 |(— + H.c. 

and Pai(0) denotes asymmetric interference part of the density operator in which one of the 
modes is in a diagonal coherent states basis while the other is in off-diagonal basis or vice 
versa and reads: 

Pai(0) = lexp(-i</>i) |ai)|a 2 )(a 2 |(-ai| + K)| - a 2 )(-a 2 |(-a;i| 



(11) 



H.c. 



+ exp(-i0 2 ) |ai)|a 2 )(-a 2 |(ai| + | - a\)\a 2 ) {-a 2 \{-ai\ 

The quantum properties of the system can be traced via the evolution of the normally 
ordered characteristic function, which for the two-mode case has the form 

CV(G,C 2 ,t) = Tr |p(0)nexp[Oi](t)]exp[-C;^(t)] 

where, for the system under consideration, Aj(t) are given in (jlj) and p(0) is the initial density 
operator of the system ([7]). It is obvious that the characteristic function here includes 16 ele- 



(12) 



ments 



341 ] . We give only the calculations related to the element exp(z0i+z</) 2 )|o;i)|o; 2 )(Q! 2 | (a^J 
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of the density matrix. On substituting this element together with (TjJ into (1T2|) and calcu- 
lating the expectation values, we arrive at 



H(iX2,t) = exp 



2 



x 



(13) 

exp 1 Gc 2 £(t) + cra^* it) - [t iQi^Mt)} + e icMt) - c*«^)] 

the explicit forms of the quantities Bjj^{t) , D(t) , D*(t) and are given in the appendix 
B. It is worth mentioning that throughout the calculation of (Tl3|) , the standard commutator 
of the Langevin forces (EI) and the usual Wigner-Weisskopf approximation have been used, 
which leads to a replacement of the coupling constants by the cavity decay rates jj > 0. 



More details about the calculation of (Tl3|) can be found in [30|, l31( . On the other hand, the 
characteristic function of the single-mode case can be obtained from that of the two-mode 
case by simply setting the parameter related to the absent mode equal to zero, e.g. the 
characteristic function of the signal mode (first-mode) can be obtained by setting ( 2 = in 
(TI2]) and ([33). 

The different moments of the bosonic operator system can be determined by differentia- 
tion of the normally ordered characteristic function through the relation 

(\[Ap{t)A n /{t)) = II g ^ ^^)k^. (14) 

In order to investigate squeezing property for the compound-mode case we can define 
the two quadrature operators X(t) = \ Y^A(t) + AUt)), Y(t) = ± J2[Aj{t) - A](t)), 

where [X(t),Y(t)\ = i\ and then the uncertainty relation reads ((AX(t)) 2 ) ((AF(t)) 2 ) > \, 
where, e.g. {(AX(t)) 2 ) = ((X(t)) 2 ) — (X(t)) 2 . Therefore, we can say that the system is 
able to yield two-mode squeezing if the squeezing factor S(t) = 2((AX(t)) 2 ) — 1 < or 
Q(t) = 2((AF(t)) 2 ) - 1 < 0. Similar quantities can be defined for the single-mode case. 
We proceed to the two-mode normal generating function which is defined as 

2 

C 8°M = J^xfj J exp(-iXJl0l 2 )^(Ci,C2,t)rf 2 Cirf 2 C2. (15) 

This function may be used in studying the sum photon-number distribution and the reduced 
factorial moments for compound modes. For the general term ffl3|) . relation f fl5|) can be 



s 



calculated to obtain 



where 



i 



(1+AA+)(1+AA_) 



exp 



A+\ A_A 
1+AA+ ' 1+AA_ 



x exp[i0i + i<p' 2 — \ XXK'I 2 + \ a 'j\ 2 - 2o*«j)], 



A± = i[Bur(t) + B 2M{t)] ± y[B lM {t)-B 2M (t)f + A\D(t)\< 



(16) 



±l {a[{t)a' 2 {t)D{t) + a x {t)a 2 (t)D*{t) 



A_-A4 



(17) 



- ai(t)ai(t)[52^(*) - A±] - ^(^(^[^(t) - A±]} . 



Expression (llBTl shows that each term of the generating function of the system is the two- 
fold generating function for Laguerre polynomials. Further, in each term the quantities 
A± and \± play the role of the mean numbers of coherent photons and mean numbers of 
chaotic photons, respectively. The sum photon-number distribution P(n, t) and the reduced 
factorial moments (W k (t)) for the compound- mode case can be defined by means of the 
derivative of (A,t) via the relations 

d k 



P(n,t) 



On using (TIB]), these quantities can be deduced to become 
P/(n, t) = ( 1+A _)( 1+A+ ) exp j^- + j^t- 



A=0 



(18) 



x 



exp[z0i + i02 - | E(l"il 2 + l"j| 2 - 2a j a 'j)] 

3=1 



(19) 



X Yl (n-W 



1=0 



(n-l)U\ \1+A 



A_ V 1 - 1 ( X+ ^ l 



1+A+ 



L n -j[ 



t r ^ 

a_(i+a_)J 1j Ua + (i+a + )J' 



<w?(*)> = exp^ + ifa - § £(Kf + l«;-l 2 - 2«*«;.)] 

3=1 



X 



A_ / ' V A-l h 



(20) 
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where the subscript I in ([19]) and (!20|) means that these quantities have been calculated for 
the general term f JT3|) . Also in f fl9|) and (120]) L n (.) are the Laguerre polynomials of order n 
having the form 

Mx) = ± 7 m^. (21) 



m=0 



(n — m)!(m!) s 



We should mention that the phase terms, e.g. the term in the second line of (119]) . exist only 
in the off-diagonal elements indicating that these elements are suppressed faster than the 
diagonal elements. In other words, the greater the "distance" between the components of 
the cat states, the more rapidly the off-diagonal elements are damped. 

On the other hand, the corresponding quantities for the single-mode case can be obtained 
with the help of the single-mode normally ordered characteristic function and applying the 
same procedures as before. To be more specific, a general term J'(Ci, t) of the first-mode case 
can be obtained from ffTB"]) by setting C2 = and consequently the corresponding photon- 
number distribution and the reduced factorial moments are 

2 

2 



Pi'(ni,t) = exp 



id* + j(h> - "^KM 



3=1 



(22) 



x 



"1 



1 



,L ni [ 



!+%(*) niL B 1A r(t)(l+Buv (*))•" 



(W#(t)) =exp 



(23) 



Finally, the quantum properties of the system can be visualized well by analyzing its 
Wigner function (W). Actually, this function is sensitive to the interference in phase space 
and therefore it is helpful to be examined. The single-mode Wigner function is defined as 

W(z, t) = vr- 2 J d 2 Ci exp(-i |Ci | 2 )C^(Ci, t) exp« - Ci**), (24) 

where CV(Cij^) is the single-mode normally ordered characteristic function. The Wigner 
function for the single-mode general term is 

2 

2 



W r (z,t) 



7r[l+2fl 1Ar (t)] 



exp[i0i +i(p' 2 -\ Y.(\ a j\ 2 + \ a 'j\ 2 - 2a j a 'j)} 
3=1 



(25) 



x exp 



? W(«M[aiW-«] 

l+2Buv(t) 



}■ 
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In this expression By^(t) is always positive and represents the sum of quantum noises in 
the field, related to quantum fluctuations in the interaction and the mean value of reservoirs 
oscillators (the number of photons or phonons transferred from the reservoirs to the quantum 
system). This shows that the width of this part of W function becomes broader when the 
coupling with environment is considered. 

In the following sections we analyze the behaviour of the system under discussion on 
the basis of the results of the present section. We investigate the properties of both the 
single-mode and compound-mode cases. 



III. RESULTS FOR THE SINGLE-MODE CASE 



In our analysis for the dissipative case we consider symmetrical losses that is 7j = 7 
and (hjd) — n for j — 1, 2. We analyze two cases: underdamped case when 2g > 7 and 
overdamped case when 2g < 7. The origin of these conditions can easily be recognized by 
careful examination of the time dependent coefficients in appendix A. For example, it holds 
that 

Mt) = \ {exp[(s - l)t] + eM~(9 + ■ (26) 

In the following we discuss the behaviour of W function, photon-number distribution, 
squeezing and reduced factorial moments for the single-mode case. 



A. Wigner function 

Investigation of the single-mode W function is important since it is informative and 
sensitive to the interference in phase space, and also it gives a prediction to the possible 
occurrence of the nonclassical effects in the system. Furthermore, this function can be 



35). 



obtained experimentally using the optical homodyne tomography 

In general, the W functions of ECS, OCS and YSS (at t = 0) are consisting of two Gaus- 
sian bells, resulting from statistical mixture of individual composite states and interference 
fringes in between (signature of the nonclassical effects) originating from the superposition 
between different components of the states. Nevertheless, the locations of the extreme values 
of the interference fringes of ECS, OCS and YCS are quite different. There are several pa- 



pers, e 



•g., HQ 



18 



361 ] that have been devoted to these fringes making them less or more 
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pronounced by allowing the cat states to evolve in different quantum optical systems. For 
example, the decay rate of the interference fringes of the W function of cat states of a mode 
which is coupled to a phase- sensitive reservoir (squeezed reservoirs) 12L 34j p. 114, can be 
smaller or larger than that for ordinary thermal heat bath reservoir Uj. Furthermore, in 
these systems the fringes can be washed out completely for a specific choice of the interaction 
parameters. In the language of the density matrix this means that the contribution of the 
off-diagonal elements is small compared with the contribution of the diagonal elements but 
it is not necessarily zero. Of course, such diagonalization produces a statistical mixture of 
coherent states, which are close to classical wavepackets. Now if we consider the W function 
of the signal mode for the system under discussion, the situation becomes more complicated 



than before 
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18] owing to the entanglement between the signal and idler modes. The 



structure of the density matrix (jSJ) carries more information, e.g., when the signal and idler 
modes are initially prepared in a distinguishable macroscopic cat state and after switching 
on the interaction, the initial two Gaussian bells of the W function of the signal mode can 
be transformed, in principle, into four-fold form as indicated in ([9]), and the initial inter- 
ference fringes of the cat will be dramatically changed during the interaction as a result of 
the competition between the different components of Pai(O) and psi(O). Moreover, during 
evolution, the rate of movement of the centres of peaks of the W function is rather different 
and this leads to irregularity in its behaviour. So that the evolution of the cat states in the 
present interaction yields different types of multicomponent cat states {37I ■ 
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(a) (b) 

FIG. 1: The W function of the signal mode when the two modes (signal 
in ECS for t = 0.55, g = l,ipi = ip2 = 0,7 = h = and = § and for: a) 
(M,|a 2 |) = (3,2); c) (|ai|,|a 2 |) = (2,3). 



(c) 

and idler) 
(|ai|, \a 2 \) 



are initially 
= (2,2); b) 
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All these facts can be seen in Figs, la-c, where we have displayed the W function for 
the signal mode (excluding losses) for the given values of the parameters. In these figures 
(and throughout this paper) we have taken (Xj = \aj \ exp(iipj), where ipj is the phase of the 
initial amplitude of the jth mode and also we have defined z in (12"31) as z = x + iy where the 
quadrature variances ((AX(t)) 2 ) and ((A7(t)) 2 ) of the signal mode are associated with the 
real part x and imaginary part y, respectively. It is worth mentioning that the W function 
curves contain most of their information within one plane, in particular, for constant values 
of y. For this reason and for simplicity we have plotted a W function in one plane considering 
y ~ —0.25, where the W function illustrates more pronounced nonclassical effects (e.g., by 
means of negative values) for the case |ai| > |ct2| and also it includes a representative 
information for the other cases. Fig. la shows that the initial nonclassical negative values of 
the W function have been smeared out and the system collapses to four-component statistical 
mixture state. Further, such behaviour reveals that decoherence can be established via the 
entanglement between different modes in the parametric processes. In these cases the single- 
mode behaviour undergoes amplification resulting from the spontaneous pump photon decay 
. It should be borne in mind that in the present case the system is completely isolated and 



then such type of decoherence can be called nondissipative decoherence [39j]. On the other 
hand, the nonclassical negative values can be recovered by controlling the "distance" between 
the initial cat states of the signal and idler modes. To be more specific, when |ai| > ja^l these 
negative values can be realized, however, when |a 2 | > | on | they will disappear as indicated 
in Figs, lb and c, respectively. In Fig. lc a three-peak structure is dominant. Comparison 
of Figs. la,b and c is instructive. So, excluding the influence of the environment the system 
can decohere and recohere by adjusting the initial distance between the components of the 
cat of the input modes. However, for the dissipative case the decoherence process is related 
not only to the amplification of the pumping field but also to damping of radiation caused 
by the flux of coherent energy from radiation to the reservoirs and noise from reservoirs to 
the radiation. Under these circumstances the decoherence can be achieved in a time shorter 
than that for the undamped case. We displayed the W function for the overdamped and 
underdamped cases in Figs. 2a and b, respectively, for the same situation as in Fig. lb. 
These figures show overall distortion due to the dissipative nature of dynamics. It is clear 
that the origin of the main contribution is in the diagonal elements of the density matrix 
and then the negative values of the W function of the undamped case (Fig. lb) are washed 
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(a) (b) 

FIG. 2: The W function of the signal mode for the same situation as in Fig. lb except n = 1 and 
a) overdamped case 7 = 2g + 3; b) underdamped case 7 = 2g — 1. 

out. From Fig. 2a one can observe that the W function exhibits the well-known shape for 
the thermal light, i.e. the system exhibits Bose-Einstein statistics (super-classical light). 
However, from Fig. 2b (underdamped case) we can observe that the two-peak structure is 
dominant. This situation is similar to that of the single harmonic oscillator interacting with 
the thermal bath in which a double Gaussian structure with missing oscillatory behaviour 
occurs [ll, 



12 



361 ] . The explanation of the behaviour of W function in Fig. 2 can be 



understood by analyzing the behaviour of the W function for the single-mode general term 
(|25l) . Let us restrict ourselves to the overdamped case for which the values of the parameters 
dj(t) are exponentially decaying (this can be easily checked) when the interaction is going 
on and accordingly the centres of the peaks move toward the origin and eventually the peaks 
merge with each other. The opposite situation occurs for the underdamped case. It is worth 
mentioning that in these cases (overdamped and underdamped cases) the contribution of 
the off-diagonal elements corresponding to Pai(O) is suppressed faster than that of psi(O). 
This point will be discussed in the sum photon-number distribution in the following section. 

In general, if the interaction between different components in the system occurs (re- 
gardless whether the interaction with environment is considered or not), the decoherence 
gradually increases and thus the system evolves into a mixture (Ti pfit) < 1, Piif) is the 
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reduced density matrix of the signal mode), i.e. there occurs a destruction of the inherent 
nonclassical effects. The rate of destruction is sensitive to the nature of both the reservoir 
and parametric processes. Furthermore, we should mention that the P function possesses 
somewhat similar behaviour as the W function in such a case, i.e. it can take on negative 
values in certain regions for certain values of the parameters and therefore it cannot behave 
like a classical probability distribution function. Such an effect is independent of the type of 
cat states, which are initially used in the interaction. Finally, we should stress that it is dif- 
ficult to analyze the behaviour of the off-diagonal elements to obtain concrete information, 
so we have basically concentrated on the computer simulation. 

B. Photon-number distribution 

The concept of photon is an integral part of the modern description of light and the 
discrete nature of light can be demonstrated by a photon detector based on the photoelectric 
effect. 

On the other hand, one of most interesting nonclassical effects emerging from the su- 
perposition principle is the oscillatory behaviour of the photon-number distribution. In 
general, such behaviour is closely related to the behaviour of the W function, however, this 
is necessary but not sufficient condition. For example, the photon-number distributions 
of ECS, OCS and YSS are completely different; whereas those of ECS and OCS exhibit 
pairwise oscillations in phase space (even number of photons can be observed for ECS and 
odd numbers for OCS), the distribution of YSS is a Poissonian even though the behaviour 
of the W function for these states is qualitatively similar. In the interaction under discus- 
sion for undamped case the oscillatory behaviour in the photon-number distribution can 
be established even if the initial cat states exhibit Poissonian statistics. This is of course 
based on the values of the interaction parameters. The origin of such behaviour is in the 
interference in phase space where the photon-number distribution of input coherent light is 
always displaying a single-peak structure, which is broader than the corresponding Poisson 
distribution with the same mean photon number. 

We start our discussion by investigating the behaviour of the undamped case. We con- 
sider here the photon-number distribution P(rii) of the signal mode when the signal and 
idler modes are initially prepared in the YSS. We have seen that this quantity can ex- 
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FIG. 3: Single-mode photon-number distribution of the signal mode when the signal and idler 
modes are prepared initially in YSS; the values of the parameters are the same as in Fig. lb and 
ifi = f (solid curve), — | (dashed curve). 

hibit oscillatory behaviour after switching on of the interaction by a suitable time provided 
that |o(i| > | «2 1 , as indicated in Fig. 3 for shown values of the parameters. In this figure 
ip — (j) — ipi — ip2, where ipj is the phase of the initial jth mode (signal or idler) and is 
the phase of the pump as before. In this case ip represents the phase mismatch. Further- 
more, one can see that when tp changes from n/2 to — 7r/2 the parity of oscillations changes 
(compare dashed and solid curves). The reason for taking if> = ±tt/2 can be found in 27 ]. 
The origin of the oscillatory behaviour in the photon- number distribution P{n\) here is the 
competition between the contributions of /5m (0) and of psi(O), as we will see in the sum 
photon-number distribution. Further, it is worth mentioning that the behaviour of P{n\) in 
ig. 3 is similar to that of initially cat state of a mode coupled to a phase-sensitive reservoir 



12J, [34J p. 114, however, the source of the oscillations in jl2j is the phase information 



included in the reservoir, which can be transferred to the field. In Fig. 4 we have displayed 
-P(ni) for the same situation as given by the solid curve in Fig. 3, but for |ai| < |ct2 1 (solid 
curve); further the damped cases are considered: underdamped case (short-dashed curve) 
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FIG. 4: Single-mode photon-number distribution of the signal mode when the signal and idler 
modes are prepared initially in YSS, the values of the parameters are the same as in Fig. lc with 
V> = f for 7 = n = (undamped case-solid curve); n = 1,7 = 2g — 1 (underdamped case-short- 
dashed curve) ;n=l,7 = 2<7 + l (overdamped case- long-dashed curve) . 

and overdamped case (long-dashed curve). Comparing the solid curves in Figs. 3 and 4, we 
see from Fig. 4 that the oscillations in P{n\) are smoothed out. Further, the comparison of 
various curves in Fig. 4 shows that the behaviour of P(ni) for underdamped and undamped 
cases is similar in the sense that they include smooth oscillations. These smooth oscillations 
are completely washed out for the overdamped case, as is expected. Actually, for \a 2 \ > \a±\ 
coherence is lost and the main contribution is related to the energy of the field mode (di- 
agonal terms). As we can see the behaviour of the photon-number distribution is in a good 
agreement with that of W function. Finally, it has been verified that the behaviour of P(rii) 
for the damped case when \a±\ > \a 2 \ and |aii| < |a 2 | is quite similar. This is connected with 
the fact that the oscillatory behaviour in the photon-number distribution is highly sensitive 
to the dissipation dynamics. 
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C. Single- mode squeezing and reduced factorial moments 

Squeezing is one of the most important phenomena in quantum optics because of its 
applications in various areas, e.g. in optics communication, quantum information theory, 
etc. [40]. Squeezed light can be measured by a homodyne detection in which the signal is 
superimposed on a strong coherent beam of the local oscillator. Furthermore, quite recently 
it has been shown experimentally that there is an evidence of squeezed light in the biological 
systems (41] . So an analysis of squeezing phenomenon in quantum optical systems is an 
important topic. 

We start our investigation by determining the behaviour of the single-mode squeezing. 
Generally, in the system under consideration the output mode loses its initial squeezing 



feature during the interaction 16| as a result of an amplification process and the interaction 
with the environment, which accelerates the loss in quantum fluctuations. Here we analyze 
the influence of different types of cat states on the behaviour of quadrature squeezing. To 
do so, considering aj to be real, we write down the quadrature variance of the ^-component 
of(which is expected to yield squeezing) for the signal mode and various initial input cat 
states as follows: 

Q ee (t) = X - {B lN {t) + a\fl{t){t&n\ial - 1) + a 2 2 \f 2 (t)\ 2 [t a nha 2 2 + cos(20)]} , (27) 
Q oe {t) = \ {B lN {t) + a$fi(t)(cothal - 1) + a 2 |/ 2 (t)| 2 [tanha 2 + cos(20)]} , (28) 

Qey(t) = \ {Bw(t) + a 2 f!(t)(t&nhal - 1) 

(29) 

+ a 2 2 \f 2 (t)\ 2 [l + cos(20) - 2exp(-4a|) sin 2 0]} . 

In these expressions the subscripts ee, oe and ey stand for the initial (signal, idler) modes 
which are in (ECS, ECS), (OCS, ECS) and (ECS, YSS), respectively. Furthermore, one 
can note that the significant value of squeezing can be obtained when the pump phase is 
cf) = ±7r/2. Also it can be mentioned that the correlation between signal and idler modes 
does not occur since cross terms such as a\a 2 are absent. Actually, this is connected with the 
two-photon nature of ECS and OCS where (a™(0)) = 0, j = 1, 2, when m is an odd integer. 
It is worth mentioning that the origin of losses caused by reservoirs in expressions ( 127]) -( 129]) 
is the mean photon number. Further, it should be reminded that ECS and YCS can exhibit 
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normal squeezing, which is more pronounced for ECS, however, OCS are unsqueezed states 



Further, for the undamped case with = ir/2, expression (1281) can be rewritten as 



The last term in this expression is non-negative and can become zero by appropriately 
choosing the value of a\. In order to get squeezing in the signal mode, expression (|30|) 
must be negative and this depends on the behaviour of the function f(x) = x(tanhx — 1). 
In other words, squeezing can be established if 1 + f(x) < 0, however, —0.3 ^ f(x) < 0, 
and consequently squeezing cannot be obtained. More specifically, to obtain single-mode 
squeezing from this device, the mode under consideration should be prepared initially in a 
squeezed-cat state regardless of the type of the cat state in the free port (if it is squeezed or 
not). Furthermore, one can easily estimate how long the single- mode squeezing of the initial 
light can survive if this light is imposed at the input of the parametric amplifier. Restricting 
ourselves to the input ECS and using the fact that squeezing is surviving if Q ee (t) < 0, the 
time range over which such a situation occurs is 



where f(x) has the same expression as before. For Q oe {t) (i.e. (128]) ) the procedure shows 
that the required time is a complex number and this agrees with the previous remark. 

From (I2T|) and (I29j) one can easily verify that the amount of squeezing available in Q ee {t) 
is much larger than that in Q ey {t) provided that ai is finite; nevertheless, when the value 
of «2 is zero or large enough, both quantities are typical. This implies that the larger the 
degree of squeezing in the free port is, the more squeezing will be available in the mode 
under consideration. 

Fig. 5 has been displayed to illustrate the behaviour of the Q-squeezing factor for the 
undamped case when the signal mode is prepared initially in ECS and the idler mode in 
ECS (solid curve), YSS (short-dashed curve) and OCS (long-dashed curve) in dependence 
on q.2 for the given values of the parameters. Furthermore, the triangle- and star-centred 
curves represent the underdamped case with a zero-temperature (n = 0) and a nonzero- 
temperature heat baths (n ^ 0), respectively, when the two modes are in ECS. The choice 
of a \ = 0.7 is related to the fact that the initial ECS gives maximum squeezing at this 



15] 





(31) 
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FIG. 5: The single-mode squeezing factor Q(t) for the signal mode being initially in ECS and the 
idler mode being in ECS (solid curve), YSS (short-dashed curve) and in OCS (long-dashed curve) 
for t = 0.2, g = l,af = 0.7, ipi =^2 = 0,7 = n = and (j) = |. Triangle- and star-centred 
curves are given for the squeezing factors of a zero-temperature heat bath (n = 0) and a nonzero- 
temperature heat bath (n = 0.1), respectively, for the underdamped case with 7 = 2g — 1.6 and 
for the same situation as represented by the solid curve. 

value. Now the analytical facts discussed above are remarkable in Fig. 5. Further, from this 
figure we can also see that for the OCS-idler-mode squeezing is minimum when a| is close 
to zero (indeed there is a singularity at a 2 = related to the nature of the OCS), increasing 
monotonically as a 2 increases, then it stagnates at large values of a 2 yielding its maximum 
value. Actually, the behaviour of this case is quite different from the behaviour of the ECS- 
idler- and YSS-idler-mode cases (compare long-dashed curve with solid and short-dashed 
curves in Fig. 5), where the maximum value of the former is the minimum value for the 
latter. This means that for large values of a| the signal mode can produce the same value 
of squeezing regardless of the type of the input cat states in the free port. Also this figure 
shows how one can control the single-mode squeezing relying on the type of initial cat states 
in the free port. 
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We now turn our attention to the damped case. As we mentioned earlier By^(t) in 
expressions fl27l) - fl29|) is always positive and therefore the coupling of the system with the 
environment degrades the amount of the single-mode squeezing. We have given two examples 
in Fig. 5 for the underdamped case (triangle- and star-centred curves). Comparing these 
curves with the solid one in the same figure, the conclusion becomes clear. Moreover, from 
the comparison between triangle-centred curve and star-centred curve we can conclude that 
in nonzero-temperature heat bath case the quantum coherence (i.e. nonclassical effects) 
is lost much faster than in zero-temperature heat bath case On the other hand, the 
occurrence of squeezing in the damped case has been verified also in the behaviour of the 
W function where for the same values of the parameters we obtained that it always exhibits 
noise-ellipse forms of cuts together with single-peak or two-peak structure according to «2 
is small or large, respectively. It is worth mentioning that for the cat states interacting with 
heat bath the authors of [5| claimed that squeezing is much more robust with respect to 
damping than the oscillations in the photon-number distribution or the interference in phase 
space described by W function. They obtained this conclusion by analyzing these quantities 
graphically. Nevertheless, in their analysis, squeezing has been obtained for a specific region 
of a (in particular when a is small) whereas the behaviour of the photon-number distribution 
and the W function has been analyzed in a different region (when a = 2, where squeezing 
does not occur, see Fig. 8a in jlj]). This leads to the incorrect conclusion. We have examined 



Dution and the W function for these cases 
sj], however, we have taken a = 1 (where 



the behaviour of both the photon-number distri 
using the same values of the parameters as in 
squeezing is noticeable) and found that the behaviour of these quantities reflects equally the 
properties of squeezed light. For more details reader can consult 42]. 

Finally, for the single-mode reduced factorial moments, which can be measured by a set of 
photodetectors, we have found that the behaviour of these quantities is in a good agreement 
with the behaviour of the single-mode squeezing. 



IV. RESULTS FOR THE COMPOUND-MODE CASE 

In this section we demonstrate the two-mode properties for the system under discussion by 
determining the sum photon-number distribution, two-mode squeezing and reduced factorial 
moments. 
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A. Sum phot on- number distribution 

The oscillations in the joint photon-number distribution might be observed in experiments 
that generate two-mode squeezed light in which the two modes may be distinguished by 
frequency or by propagation direction. Having separated the two modes, as was done by 
using a polarizing beamsplitter 43| , one can send them directly onto separate photocounters 



and then build up the joint photon-number distribution from the photocount statistics of a 
sequence of pulses 1441 1 . 

n 

It has been shown [27[ that the sum photon-number distribution P(n) of the nonde- 
generate parametric amplifier, when the modes are initially prepared in coherent state and 
under certain conditions, displays two regimes, which exhibit either single-peak structure 
or oscillatory behaviour. The single-peak structure can be narrower (or broader) than that 
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FIG. 6: Sum photon-number distribution when the signal and idler modes are prepared initially 
in ECS with |ai| = 3, |ck2 | = 2,g = 10 4 , i = 3 x W~ 4 ,ip = ^. On the right corner of the figure the 
scheme shows the parity of the photon-number sum. 



of the corresponding Poissonian distribution showing nonclassical (or classical) effects. The 
striking feature is that the P[n) can exhibit, for certain choice of the interaction parame- 
ters, collapses and revivals in the photon-number domain somewhat similar to those known 
in the JCM. For the present system the situation is rather complicated regarding to the 
evolution of the density matrix (JSJ) where the interference in phase space is established. 



22 



Excluding dissipation and using tjj = ±7r/2 as in [27j], we have seen generally that P(n) 
exhibits always oscillatory behaviour irrespective of which type of cat states has been con- 
sidered initially. Moreover, P(n) can yield pairwise oscillations identifying that even or odd 
photons are being observed. Figs. 6 and 7 have been plotted to show such a phenomenon 
for given values of the parameters. From Fig. 6 we observe that the long scale oscillations 



in the 

states 



aehaviour of P(n) with P(2n + 1) = are somewhat similar to those of squeezed 



45]. Actually for squeezed states the pairwise oscillations are explained as a di- 



rect consequence of the quadratic, or two-photon nature of the squeeze operator S(r), i.e. 



S(r) = exp[|(a 2 — a) 2 )\ 45 1. However, here the origin of the pairwise oscillations in P( 



n 



is in the competition between the processes described by three parts of the density matrix 
of the field. Figures 7a-c give insight into this point, i.e. they show the manner in which 
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FIG. 7: Sum photon- number distributions corresponding to the three parts of the density matrix 
for the same values of the parameters as in Fig. 6 and for: a) P(n) resulting from /5m(0); b) P(n) 
resulting from /Ssi(0); and c) P(n) resulting from pai(0). 



the photon-number distributions of the three parts of the density matrix ([8]) compete. More 
precisely, Figs. 7a,b and c represent the sum photon-number distributions Pm(u), Psi{n) 
and Pai(^) associated with p M (0),p SI (0) and Pai(O), correspondingly. Fig. 7a shows revival- 
collapse pattern, which is resulting from the statistical mixture part (Q. Furthermore, from 
this figure one can observe that when the number of photons increases, the amplitude of 
the revivals diminishes, but the revival periods extend. It is worth mentioning that such 
a quantum collapse-revival phenomenon has been seen for the photon-number distribution 
of single-mode [46 1 and two-mode [42 1 squeezed coherent states with complex squeeze and 
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displacement parameters. Associations with the quantum phases of the modes may be con- 



sidered 
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481 ] . 34J p. 281. From Fig. 7b it is clear that Psi{n) oscillates between negative 



and positive values. Comparison of Fig. 7a, Fig. 7b and Fig. 7c shows that the values of 
Pai(ti) are approximately negligible compared with those of -Pm(^) and Psi{n). Now if we 
turn our attention back to Fig. 6 we can recognize that the competition between Pm(^) and 
Psi{n) leads to the destruction of collapse-revival phenomenon (which appeared in Fig. 7a), 
however, the distribution "evolves" in an interesting manner similarly as that for squeezed 
states. In fact, the situation here is in contrast with that of the evolution of cat states in 
the JCM where the interference in phase space makes the revival- collapse phenomenon in 
the atomic inversion more pronounced when the revival time equals the half of that for the 
standard JCM with initial coherent light [l?J] . The reason is that for the latter case the 
distribution of the spectral components of the atomic inversion is one dimensional distribu- 
tion since only one mode is involved, whereas here we have two-dimensional superimposed 
distributions. Further, we have obtained that the oscillations in P(n) increase if either g or 
t or both are increased. This fact is clear if we look at the problem as evolution of cat states 
under the action of two-mode squeeze operator, where the squeezing parameter in this case 



is r = gt and the oscillations become more pronounced for the large values of r [44! . 149 1. 
In conclusion, for the system under consideration the origin of the oscillations in the sum 
photon-number distribution is two-fold: (i) The interference in phase space, (ii) The strong 
coupling between the signal and idler modes of the system in the course of the interaction 
time. 

On the other hand, we have found by an explicit algebraic calculation for all quanti- 
ties studied in this paper excepting quadrature squeezing that the contributions associated 
with elements of psi(O) and Pai(O) involve cos(0i ± 4> 2 ) and cos0i i2 , respectively. This fact 
together with the information included in Figs. 7 show that P(n) can collapse and yield 
decoherence in relation to distribution for statistical mixture part using specific types of cat 
states initially. More illustratively, preparing one of the modes initially in ECS (or OCS) 
and the other in YSS or vice versa, the P(n) evolves as described by Fig. 7a (for the same 
values of the parameters). This shows how one can decohere the system apart from the 
amplification nature of the system and without coupling it to the environment. We call such 
a type of decoherence a phase decoherence. Actually, for some quantities the contribution 
of off-diagonal elements of the density matrix is responsible for the nonclassical effects, e.g. 
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as we will see below in the compound- mode reduced factorial moments. Then we can make 
a good estimation of these phenomena using such a property. The final remark is that the 
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FIG. 8: Sum photon-number distribution for a) undamped case with the values of the parameters as 
those of the solid curve in Fig. 3; b) damped case with (g, n, 7) = (0.5, 0.5, 2g — 0.9) (underdamped 
case-solid curve) and (0.5, 0.5, 2g + 0.1) (overdamped case-dashed curve). 

nonclassical effects for the compound-mode case are much richer than those for the single- 
mode case. This results from the strong quantum correlation between the signal and idler 
modes, which manifests itself as the summation in expressions (jl~9j) and (|20|) . To make this 
point clear, we give Fig. 8a for the P{n) and for the same situation as that of the solid curve 
in Fig. 3. The comparison between these two figures is instructive. Further, Fig. 8b includes 
information on the damped case where one can observe smooth oscillations for the under- 
damped case (solid curve) and two-peak structure for the overdamped case (dashed curve). 
Specifically, for the overdamped case the off-diagonal elements are completely suppressed. 



B. Two-mode squeezing and reduced factorial moments 

Firstly, we will discuss two-mode squeezing, which can be measured by heterodyne de- 
tection where squeezing is carried jointly by two modes of different frequencies and the local 
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oscillator has a frequency midway in between [50|. We proceed taking into account that 
the nondegenerate parametric amplifier is well described by the two-mode squeeze operator. 
Actually, two-mode squeezing phenomenon is essentially connected with the states upon 
which such operator acts. In other words, squeezing may not exist for some specific states 
even if they include two-mode squeeze operator in their structures. To show this we give 
the form of the two-mode squeezing factor Q when the signal and idler modes are initially 
prepared in cat states, excluding losses, considering aj,j = 1,2, are real and = 0. In this 
case the squeezing factor Q takes the form 

<?(f) = ~[Qi(f) + <&(*)], (32) 

where Qj(t),j = 1, 2, represent the corresponding single-mode squeezing factors of the signal 
and idler modes, respectively. From expression (1321) it is clear that the correlation between 
the signal and idler modes does not exist since {YiY 2 ) = (^i)(^2), where (Yj),j = 1,2, are 
the expectation values of the F-quadrature of the signal and idler modes. The existence 
of correlation between these two modes is important to obtain squeezing: in the compound 



modes even if the individual modes are not themselves squeezed [49] . Expression (1321 shows 
that interference in phase space can destroy squeezing and also that the rates of degradation 
of the two-mode and single-mode squeezing are on the same level for symmetrical losses. 
Also it is obvious that the two-mode squeezing can be realized if at least one of the two 
modes (signal or idler) can exhibit single-mode squeezing. This is necessary but not sufficient 
condition. Further, the maximum squeezing may be produced when both the signal and idler 
modes exhibit maximum single-mode squeezing. This reflects the importance of the choice 
of the type of the initial cat states. From the discussion in section 3 it is obvious that when 
the two modes are initially in OCS, the interaction cannot generate two-mode squeezing. 

We now discuss the case that one of the modes is initially squeezed (ECS) and the other 
is unsqueezed (OCS) for the same situation as that in (|32|) . The required time to obtain 
squeezing in the compound-mode case (Q(t) < 0) is 
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This inequality gives relation between the squeeze time and the "distances" between the 
states, which are forming the cats. So that the length of the device can be adjusted to 
obtain squeezed light. An analysis to the right-hand side of f )33|) gives that the maximum 
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FIG. 9: The compound two-mode squeezing against = 1,2, for undamped case when the 

signal and idler modes are initially prepared in ECS for t = 0.2, \aj\ = 0.7, g = 1,0 = §: a) 
squeezing factor S; b) squeezing factor Q. 

scaled time (effective time) for squeezing is r = gt ~ 0.1769 at (|ai|, |ok 2 |) = (0.6,2.4). It is 
worth mentioning at this point that the F-component quadrature squeezing of ECS takes 
on its maximum, whereas that of OCS vanishes. 

On the other hand, when the phase information is included (i.e. ipj ^ 0, ^ 0) the 
situation is improved in such a way that squeezing can be detected in both quadratures 
provided that one of the modes is initially prepared in a squeezed-cat state, i.e. in ECS or 
YSS. The role of the phase is analyzed in Figs. 9a and b where we have plotted the squeezing 
factors corresponding to X-component and F-component, respectively, when the signal and 
idler modes are initially prepared in ECS. From these figures we can see that maximum 
squeezing occurs in X-component and locates at (1^1,1^2) = ( m i' n ', m 2' K ), m j = 1/2,3/2. The 
opposite situation can be observed for the F-component where squeezing also exists but 
it is less pronounced and locates at (1^1,1^2) = (mi7r, m 2 7r), = 0,1,2, i.e. the extreme 
values are exchanged. Such a behaviour in the two quadratures indicates that the uncertainty 
relation holds. Further, for = —it/ 2 we obtained similar behaviour, however, the maximum 
squeezing values would be available in the F-component. Further similar behaviour can 
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FIG. 10: The compound-mode reduced factorial moments [K c = (W k ) / (W) k — 1) against |ai| 
for k = 5,tf) = n/2,t = 0.2, |ck2 | = 0.25,5 = 0.5, and for the undamped case when the signal and 
idler modes are prepared initially in (OCS,OCS) (solid curve) and (OCS,ECS) (long-dashed curve). 
The damped case is considered for the same situation as the solid curve, but for the underdamped 
case (short-dashed curve) (7, n) = (2g — 0.6, 0.5) for the overdamped case (star-centred curve) 
(7,n) = (2g + 0.1,0.5). The straightline shows the antibunching bound. 



be seen if the other types of cat states are used. It is remarkable that the behaviour of 

n 

squeezing in the present interaction is quite different from that for both ECS [3[ and two- 



mode squeezed coherent states 



49| where for both squeezing exists only in one of the two 



quadrature components; perfect squeezing (i.e. 100% squeezing) occurs only for the two- 
mode squeezed coherent states. In conclusion the squeezing phenomenon discussed here is 
a manifestation of co-operative effects among the phases of the different components of the 
system. Furthermore, by controlling the phases in this device, squeezing can be amplified 
and switched between the two quadratures. 

We would like to conclude this subsection by shedding the light on the behaviour of 
the reduced factorial moments for the compound-mode case. In this case the antibunched 
light can be measured if both the signal and idler modes are simultaneously detected by 
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means of two photodetectors and then their outputs are correlated. Generally, reduced 
factorial moments can evolve to produce nonclassical negative values based on the types of 
initial cat states. Further, we have checked that the contribution of pu(0) cannot produce 
antibunched light independently. These two facts reflect the role of the phase decoherence 
property discussed above. On the other hand, similar argument as that for the single- 
mode case can be given here except the case when the two modes are in YSS where a slight 
antibunching can be observed, but only for certain values of the interaction parameters. Fig. 
10 shows the reduced factorial moments of the compound-mode case for the given values of 
the parameters. In this figure we have used the normalized factor K c (= (W k )/(W) k — 1), 
where the nonclassical effect occurs when K c < 0. From this figure, it is interesting to 
observe that the behaviour of the reduced factorial moments when the (signal, idler) modes 
are in (OCS,OCS) is very close to the behaviour of the second-order correlation function of 
OCS [5| provided that \ocj\,j = 1,2, are finite. 



V. CONCLUSIONS 



In this paper we have analyzed the properties of the dissipative parametric amplifier when 
the signal and idler modes are initially prepared in Schrodinger-cat states. Needless to say 
the present interaction as described by the Hamiltonian ([3]) and the density matrix (jSJ) is 
much more complicated than that in the simple case of a harmonic oscillator coupled with 
the heat bath ll|- 18] . After obtaining the solution of the Heisenberg-Langevin equations 
quantum statistics of interacting modes have been investigated based on the normally or- 
dered characteristic functions. The system of damping oscillators has been assumed to have 
a flat spectrum and a chaotic distribution. 

In general, there are two operations controlling the behaviour of the interaction which are 
the interference in phase space and the entanglement. Intuitively, the initial macroscopic 
cat states cannot be preserved in the system. Furthermore, for long-time interaction the 
initial nonclassical effects of the cats are degraded by the amplification dynamics inherent in 
the system and the cumulative effects of dissipation. Excluding dissipation we have shown 
that if the input to the device are squeezed (sub-Poissonian) cat states, then the output 
may be squeezed (antibunched) too provided that the interaction time and gain are finite. 
Furthermore, the device could be used to amplify squeezing in the compound-mode case. 
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On the contrary, the well-known role of the parametric amplifier as a source of perfect 
squeezing in the compound-mode case rather than in single-mode case fades out for input 
cat states. The photon-number distribution can exhibit oscillatory behaviour, which is more 
pronounced in the compound-mode case than in the single-mode case, regardless of the types 
of the initial cat states. Also we have shown that the decoherence can arise in the system 
from the decay of the pump and the phase control. 

For dissipation case we have considered in detail only two cases which are related to un- 
derdamped and overdamped regimes. In these cases the interaction tends to eliminate the 
off-diagonal elements of the density matrix and to affect the diagonal elements. It has been 
found that the dynamical subsystems (signal or idler modes) collapse to the statistical mix- 
ture state or thermal state according to whether the underdamped regime or overdamping 
regime is considered. Further, we have also shown that the single-mode squeezing at nonzero 
temperature of the environment decreases much faster than that at zero temperature, how- 
ever, this situation is generally valid for all quantities studied in this paper. In conclusion, 
the losses uniformly distributed over the device deteriorate its operational characteristics. 
Of course, such deterioration is more pronounced for the overdamped case. Finally, the 
inclusion of lossy mechanics in the system is of a great interest for accurate measurements. 
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Appendix A 



In this appendix we give the explicit forms for the time-dependent coefficients 
(fj(t), Tji(t), r^ 7 (i)) of the solution of the Heisenberg-Langevin equations of the Hamilto- 



nian (j3J) (for details see |30L l3lj' 



(7i + 72)* 



/i( f ) = "^ ex P 
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where 



e = (71 - 72) 2 + 16c/ 2 , e' 7 = -(71 + 72) + z(u/,- - <^ 7 )- 



(A5) 
(A.Q) 



Appendix B 



In this appendix we write down the explicit forms for the quantities (Bj^(t), D(t)) in the 
expression (TT3"j) (for details see 



B m (t) = £iaK7i,72, (nu), (n M },t) = ^{8c/ 2 £i 
+ [(726 - 4/(7i + 72)) £ - V~e (72(72 " 7i) + V) F] 



7172 - 4# 2 

472^ 2 (1 + (n 2d )) ^ 
2 — [(71 + I2) E - ^JeF] 



7i72 - 4g 
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D(t) = 29eX ^ ) {( l2 - ll )E 1 - v TeF 
:[(7i72-7^-8^ + 72^] 
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(5.3) 



<*i{t) = oc\h{t) + a'j;(t), a 2 (t) = a'J*(t) + a* 2 f 3 {t), 
a[(t) = a'Mt) + a* 2 f 2 (t), a' 2 (t) = a*J 2 (t) + a' 2 f 3 (t), 



(BA) 



31 



while 

E — 1 — exp 



(7l + l2)t 



cosh ( ^-t 



i?i = exp 
F = exp 
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sinh I 
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(5.56) 
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A generalization to the case when the total density operator is not factorized and the 
correlation between the system and reservoirs is present, may be excluded 51 |. 
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